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Abstract. A quasi- local energy for Einstein's general relativity is defined by the value 
of the preferred boundary term in the covariant Hamiltonian formalism. The boundary 
term depends upon a choice of reference and a time- like displacement vector field (which 
Q^' can be associated with an observer) on the boundary of the region. Here we analyze 

5^ . the spherical symmetric cases. For the obvious analytic choice of reference based on 

W)| the metric components, we find that this technique gives the same quasi-local energy 

values using several standard coordinate systems and yet can give different values in 
some other coordinate systems. For the homogeneous-isotropie cosmologies, the energy 
' can be non-positive, and one case which is actually flat space has a negative energy. As 

■ an alternative, we introduce a way to determine the choice of both the reference and 

lO , displacement by extremizing the energy. This procedure gives the same value for the 

' energy in different coordinate systems for the Schwarzsehild space, and a non-negative 

. value for the cosmological models, with zero energy for the dynamic cosmology which 

is actually Minkowski space. The timelike displacement vector comes out to be the 
dual mean curvature vector of the two-boundary. 



! 1. Introduction 



The identification of a good expression wliicli well describes tlie energy for gravitating 
systems (and really this means for all physical systems) still remains an outstanding 
issue. One consequence of the equivalence principle is that there is no well defined 
(i.e., covariant) local description of the energy-momentum for gravitating systems (see, 
e.g., the discussion in §20.4 of [1]). The modern idea is that properly energy is quasi- 
local (i.e., is associated with a closed 2-surface, for a comprehensive review see [2]). 
One approach is to regard energy as the value of the generator of dynamical changes 
with time, the Hamiltonian. Here we consider in particular the covariant Hamiltonian 
formalism [H HI O [6] . Within that approach a certain preferred Hamiltonian boundary 
term was identified [6]. The quasi-local energy is given by the value of this boundary 
term with a suitable choice of time evolution vector field on the closed 2-surface. In 
addition to the spacetime displacement vector field and the dynamical variables, this 
boundary term also depends in general on a choice of certain reference values for 
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the dynamical variables (which represent the ground state with vanishing quasi-local 
quantities). Unlike the case for other fields, the reference values for gravity theories 
based on dynamic geometry (the metric and connection) cannot be taken to have trivial 
values. Essentially this is because the natural ground state for dynamic geometry, 
Minkowski space, has a non-vanishing metric. Hence the choice of reference for such 
theories in general, and thus in particular for general relativity, necessarily requires some 
suitable way to select an appropriate Minkowski geometry at the points of the closed 
2-surface. Presently this is a quite active research topic [HIHlEj- Here, following [TO] , 
we consider this problem for the most important special case: spherically symmetric 
solutions (more specifically, Schwarzschild and the homogeneous isotropic cosmologies) 
to Einstein's gravity theory. General Relativity (GR). We consider two approaches. 
The first we name analytic — essentially one takes the limit of the physical metric 
components in some coordinate system when the physical parameters take on trivial 
values. This type of approach goes back to [U [11], and shows that the quasi- local 
energy is reference and displacement vector dependent. That naturally raises a question: 
is there a minimum (maximum) value among all these available choices? Which leads to 
the second approach: find the optimal reference via a variational principle extremizing 
the energy. Here we show that the analytic approach can give the same standard quasi- 
local energy value for several choices of the spatial coordinate system and yet will lead to 
different energy values for different time coordinates. The value of the obtained quasi- 
local energy is not necessarily non-negative. On the other hand the energy extremization 
always gives a coordinate frame independent quasi-local energy value which is, moreover, 
non-negative and vanishing only for Minkowski space. 

The outline of this work is as follows. In Section 2, we briefly introduce the 
covariant Hamiltonian approach, which leads to the Hamiltonian boundary term that 
gives the quasi-local energy. Section 3 concerns the analytic approach to choosing the 
reference: 3.1 includes the analysis of the Schwarzschild geometry in three different 
spatial coordinate systems and the Eddington-Finkelstein and Painleve-Gullstrand time 
slicings; the homogeneous-isotropic cosmology metrics are considered in Section 3.2. In 
section 4, we use the method of extremizing the quasi-local energy to determine the 
choice of reference and also the displacement vector. Section 5 is the conclusion. 



2. The covariant Hamiltonian approach 

Our approach to quasi-local energy is via the covariant Hamiltonian formalism, which 
has been described in detail in a series of works [31 HJ [5l [6]. The construction of the 
energy expression is briefly outlined here. It begins from a first order Lagrangian 4-form 
for some /c-form fleld. 

£ = d(/^Ap-A(^,p). (1) 

The variation has the form 

6C = d(5(p Ap) + 5(pA — + — A5p. (2) 

d(p dp 
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Hamilton's principle applied to the action obtained by integrating the first order 
Lagrangian over a region leads to the pair of first-order field equations: 

= ^:=-{-l)'dp-d^A, = ^-f := - dpA. (3) 

o(p op 

From infinitesimal diffeomorphism invariance (generated by a vector field N) for C 
it follows that (|2]) must become an identity under the replacement 6 — > (where 
£js! = d^N + '•Nd on the components of form fields): 

dij^C = £n£ = d(£N(/? Ap) + £^ip A- — h — A £j^p. (4) 

dip op 

From this identity it follows that the 3-form 

H{'N):= £^ipAp-L^C (5) 
satisfies the (a particular case of Noether's first theorem) identity 

- dnCN) = £^ip A — + — A £^p. (6) 

dip dip 

Hence (when the field equations are satisfied) the 3-form 'H(N) is a conserved current. 
From an expansion of its definition it can be seen to have the form T-L(N) = 
N^Hfj, + di3(N). Inserting this expansion into the identity (j6]) gives an identity with 
terms proportional to dA^^ and iV^, which should hold at any point for arbitrary 
values of these variables. From the terms proportional to dA^^ one learns that for 
any diffeomorphically invariant Lagrangian, "H^ = —L^ip A |^ + (—1)^^ A L^p (this is a 
special application of Noether's second theorem), and thus "H^ vanishes when the field 
equations are satisfied. Consequently, the value of the conserved current associated with 
a spatial region S and vector field N, 

E(N,S) := / H{^) = [ N^H^,+ <f i3(N) = / i3(N) (7) 

really depends only on the boundary term, and hence is actually quasi-local. The 
physical interpretation of this value is that for suitable displacements on the boundary 
it represents the components of the quasi-local energy-momentum. However it must be 
noted that (just like other Noether conserved currents) one can add any total differential 
to ([5]). The resulting 3-form would still be a conserved current, but it would define a 
different conserved value. By this process one can adjust the boundary expression 
B(N) and thereby the associated conserved energy-momentum to have almost any value. 
Fortunately the 3-form ([5]) has another important role which brings the freedom in the 
choice of the boundary term B, and hence the value i?(N,9S), under physical control. 
The 3-form 'H(N) is not only the conserved current associated with an infinitesimal 
diffeomorphism it is also actually the Hamiltonian density., i.e., the Hamiltonian (N, E) 
that generates the dynamical evolution of the variables along a timelike vector field N 
is just ([7]). The variation of this Hamiltonian has the form 

(5iJ(N,S)= [ {6ip A £^p+ £^ip A6p) + <f C(N), (8) 



Quasi-local energy and the choice of reference 



4 



yielding the Hamilton equations £f^ip = 6'H(N)/Sp, £j^p = —S'H(N)/Sip. The key point 
is that requiring functional differentiability of the Hamiltonian (i.e., the vanishing of the 
boundary term in the variation of the Hamiltonian), determines the boundary conditions 
built into the Hamiltonian. Hence one should thus choose the particular form of the 
Hamiltonian boundary term B that gives the desired type of boundary condition for 
the dynamical variables (e.g., Dirichlet or Neumann) which is suitable for the physical 
problem. 

The boundary term C(N) in the variation of the Hamiltonian ([8]), will not have 
vanishing limiting value at infinity with the usual field fall offs [121 [13] — unless one 
adjusts by hand the total differential (i.e., the boundary term) in the Hamiltonian. 

Investigations led to several explicit alternative boundary term expressions. In 
general these expressions require, along with the dynamical variables on the boundary 
certain non- dynamical reference values, that represent the ground state of the physical 
system. Whereas it is generally possible, and indeed appropriate, to choose trivial 
(i.e., vanishing) magnitudes for these reference values for all the other physical fields, 
that cannot be done for dynamic geometry gravity theories, simply because the ground 
state of the metric is not a vanishing value but rather the non-vanishing Minkowski 
metric. Moreover, in a general coordinate system, the Minkowski connection also has 
non-vanishing components. 

Among the possible boundary terms corresponding to various boundary conditions, 
a certain preferred boundary term for the covariant Hamiltonian for Einstein's GR was 
identified which should be suitable for most applications. It corresponds to holding the 
metric fixed on the boundary. This choice has the virtue of directly giving not only the 
ADM quantities at spatial infinity but also the Bondi energy and flux at null infinity, 
and moreover under certain conditions it will give a positive quasi-local energy. Our 
preferred boundary term for GR is given by 

BiN) = ^{AT^, A 6,r// + ^I'^iV^lAr/^,), (9) 

where = g''^DxN^' = g''^{^xN^' + T''x-yN^), AF^^ = T^'^ - is the difference 

of the connection one-form between the dynamic space and the reference space, with a 
bar referring to the reference objects, 7]^^ = le^ti^a/?^" A 'd'^, and is the orthonormal 
coframe. 

For a given dynamical region and given dynamical fields, the value of this boundary 
expression with suitably selected vector fields can be used to determine values for the 
quasi-local energy-momentum. However, to find the specific values for these quasi- 
local quantities, one still needs to explicitly select the reference configuration and the 
appropriate displacement vector field. A natural choice for the reference is one with a 
Minkowski metric. But one must determine exactly which Minkowski space should be 
used. In the following sections we explore two approaches to achieving this. 
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3. Choice of reference: the analytic approach 

The boundary expression ([9]) includes the reference values for the dynamical fields, 
but gives no restriction as to what the reference should be. In general, the reference 
space could be any four dimensional manifold endowed with a Lorentz signature metric 
tensor and a connection. Here we take the reference space to have a Minkowski metric. 
The quasi-local energy ([7]) is the value determined by the boundary integral with a 
certain time displacement. The boundary term includes the reference variables 
for (certain projected components of) the four-dimensional metric and connection on 
the two-boundary. A reasonable requirement for choosing the reference is isometric 
embedding of the two-boundary into the chosen reference space. Without additional 
conditions, however, the embedding is not unique. 

Here we consider just spherical symmetric spacetimes. This is an important 
yet relatively simple special case for finding an isometric embedding easily, and also 
it simplifies the Hamiltonian boundary expression ([9]), so that only the first term 
i-nVii'^ contributes. 

3.1. Schwarzs child geometry 

For the static, spherically symmetric Schwarzschild metric, we will consider five different 
representations related to different time and spatial coordinates: the standard {t, r} 
where r is the area coordinate, the isotropic spherical {t, R} and isotropic Cartesian 
{t,x,y,z} coordinates, and the Eddington-Finkelstein {t,r} and Painleve-Gullstrand 
{r, r} coordinates. In each case the choice of reference is obtained analytically by 
taking m = in the metric and connection coefficients of the dynamic spacetime in 
the particular coordinates chosen. In general it can be expected that the resulting 
quasi-local energy value will depend on the choice of coordinates. 

Let us first illustrate the procedure in detail using the standard area coordinate. 

3.1.1. Standard Schwarzschild. The standard form of the Schwarzschild metric is 



The radial coordinate is determined geometrically in terms of the area of a 2-sphere: 
A = Anr"^. Take the orthonormal coframe to be 



The Levi-Civita connection one-form coefficients are obtained using the torsion free 
condition, d'd'^ + T^^, A = 0. Due to the metric compatibility condition, the 
orthonormal frame connection coefficients are anti- symmetric. The independent 
connection coefficients are 




1-— dt'+ 1-— dr2 + rW + r2sin2^d(^l (10) 
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r°i = ^dt, r°2 = o, r°3 = o. (12) 



m 

Take m = in ffTTj) and f|T2|) to obtain the reference geometry components. Then the 
non-vanishing differences of the connection components in ([9]) become 



AT\ = - ^1 - d^, Ar^s = - y 1 - sin^dy.. (13) 

Note that the term Arj^j^^ vanishes in the boundary integral for this reference choice, so 
the boundary expression ([9]) reduces to just the first term. 

Another important role in the boundary expression is played by the displacement 
vector. We assume it is a future time-like vector field. Suppose N is normal to the 
two-surface (which we choose here to be the constant t, r sphere, with 62, 63 being the 
two-surface tangent vectors), then the factor Lj^rj^'^ is obtained using N = N'^eo-\- N^ei. 

t-KiVo^ = 0, LkiTJo'^ = —r smON^dif, tN?7o^ = rN^dO, 
L^r]i^ = r sin dN^dtp, L^r]i^ = -rN^dO, 



L^r]2^ = :dr - - 2m/rN^dt. 

Jl — 2m/r 

Only the purely angular components of the quasi-local boundary term will contribute 
to the integral over the 2-sphere S: 

= 2[ArS A LNVl^ + Ar^3 A LNVl'] 



= A{1 - ,Jl -2m/r)N'^rsmede Adif. (14) 

The quasi-local energy obtained from the integral over the 2-sphere boundary at constant 
t, r then comes out to be 

Es(N) = r(l - y^l - 2m/r)iV°. (15) 

One possible choice of the displacement is the time-like Killing vector of the reference 



N = (9j = y 1 — 2m/reo, yielding 

Es{dt) = ry^l -2m/r(l - ^1^2^), (16) 
with the horizon and asymptotic limits 



Esidt)r^oo ^ mi^l - —j ^ m, Es{dt)r=2m = 0. (17) 

On the other hand, for the choice of N = (this is the unit dual mean curvature vector 
of the two-surface), the quasi-local energy value comes out to be 



^s(eo) = r(l - y'l-2m/r), (1^ 
with the horizon and asymptotic limits 



m 



Es(eo)r^oo ^ m ( 1 + — ) m, Es{eo)r=2m = 2m. (19) 
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The value f llSp is the famous result first found by Brown and York [T3]. It turns out 
that several approaches yield this same value, so we will refer to it as the standard value. 
From an examination of these two results we have noticed a curious fact: 

Es^^{dt + eo)^=m. (20) 

We do not know whether this has any significance. 

The calculations for several other metric expressions of interest to us here follow 
a similar common procedure, to avoid unnecessary repetition we have done the general 
calculation in the Appendix; we briefly report in the following the specific results. 



3.1.2. Isotropic Schwarzs child The spherical isotropic Schwarzschild metric can be 
obtained from flTII]) using the coordinate transformation r = R{1 + ^)^: 

ds' = + (l + (di?' + ^'d^' + R'^^^'Od^'). (21) 

We choose the coframe 

^° = TT^^'' ^' = (i + -/2^)'di?, 

= {l + m/2RfRde, = {l+m/2RfRsmed^. (22) 

Then we worked out the corresponding Levi-Civita connection T^^, and took m = to 
get the reference values {T^ v and 'd^). The displacement vector in the normal space of 
the two-surface has the general form N = N^Gq + A^^ei. Only the angular components 
of the quasi-local boundary term contribute. 

This procedure leads to the quasi-local energy for the spherical isotropic metric. 
The value can be calculated from expression (1A.20I) in the appendix, which gives 

Esi(N) = A°Esi(eo) = m{l + m/2R)N\ (23) 

If we choose N = (9^ (the timelike Killing field of the reference), which means 
TVO - ^-^/"^R then 

Esi{dt) = m{l-m/2R), Esi(90r^™,/2 = 0, Esi(90r->oo = m. (24) 
Another choice for the evolution vector is N = Gq, which gives 

^si(eo) =m(l + m/2i?), (25) 
with the limits 

^si(eo)ij=m/2 = 2m, £^si(eo)if^^oo = m. (26) 

The value (l25l) is just the standard result ( |T8l) . after taking into account the 
transformation r = i?(l + m/2i?)^. 

Furthermore, we considered another metric form which is also isotropic, but in 
the Cartesian coordinate system = {t,x,y,z}. This is an important check on our 
techniques, since for this representation the reference connection vanishes when we take 
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m = in the dynamic connection. With x = i?sin^cos<y9, y = i^sin^siny?, z = RcosO, 
then R"^ = x'^ + y'^ + . The metric then has the form 

ds^ = -NMi"^ + ^\dx^ + dy^ + dz^), (27) 

where TV = f^^^^/ff' $ = (1 + m/2_R)^. We choose the obvious coframe: 

= Ndt, = ^dx\ x' = {x, y, z}. (28) 

Suppose that depends on {t, R) only. Then, from the calculation given in more 
detail in the appendix, we find for the Cartesian coordinate isotropic Schwarzschild 
metric exactly the same quasi-local energy result, fl25l) . as was found using spherical 
coordinates for the isotropic Schwarzschild metric. 

Here we have calculated the quasi-local energy for the Schwarzschild metric using 
analytic matching in three different coordinate representations and obtained in each 
case the standard result. This may give one some confidence in these techniques as well 
as in the standard value. On the other hand, as we shall see in the following, there are 
other coordinate systems in which this analytic technique for determining the reference 
will lead to other values. 



ds' 



3.1.3. Eddington-Finkelstein. The Eddington-Finkelstein (EF) form of the Schwarzschild 
metric, 

2m\ 2m ~, / 2m\ , , , , , 

1 d? - 2q — dtdr + 1 + — dr^ + r^d^^, 29 

r / r V r / 

(where ^ = — 1 is for incoming and ^ = +1 is for outgoing and (t, r, 9, ip) is the standard 

coordinate system of the Schwarzschild metric) follows from the time coordinate 

transformation t = t — q2mhi{-;^ — 1), which makes the outgoing (incoming) radial 

null geodesies into straight lines of slope ±1 in the t — r plane. A principal virtue of this 

form of the metric is that it is regular at the horizon, r = 2m. Rewriting this metric in 

the ADM form, 

^ - (i . ^)-' d.. (i . 5^) (d. - c^df)^.M.^ (30) 

leads to the coframe 



# = , ^ dt, = Jl + ^ ( dr - ,^^^dl\ , 
^l + 2m/r V r \ 1 + 2m/r ) 

^'^ = rde, d^=rMd(^. (31) 

The quasi-local energy is obtained by straightforward calculating the corresponding 
Levi-Civita connection and taking m = as the reference. For N = A^'^eo + A^^ei it 
works out to be 

i?Er(N) = r I 1 - , ^ \ - ,^^=N\ (32) 

+ 2m/r ) y^l + 2m/r 
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where en = \/l + 2m/rdf + <^ — , ^"^ dr and ei = , ^ dr- For the choice of the 

V ' ^ ry/l+2m/r i/l+2m/r 

reference timehke KiUing field, N = di, this expression yields the value 

EEF{di) = 2m - r (l - (1 + 2m/r)~^/2^ , (33) 

with the asymptotic and horizon limits 

EEF{di)r^oo = fn, EEF{di)r=2m = \Plm. (34) 

For the alternative choice of N = eg, the quasi-local energy obtained from f l32p is 

EEF(eo) = r (l - (1 + 2m/r)-V2) = (35) 
with the asymptotic and horizon limits 

^EF(eo)r-^oo = £'EF(eo)r=2m = ^{2 - a/2). (36) 

From these two time choices we find the curious fact again: 

^EFQ(9i + eo)) =m. (37) 

Another choice is the unit dual mean curvature vector (outside the horizon) 

Vt-L 1 / 2m 

N = , = eo - ? — ei 

a -4mVr2 V r 



the associated quasi-local energy is 



r 



EeF^ )r^oo = fn, ^Ef(N )r=2m ^ oo. (39) 

3.1.4- Painleve-Gullstrand. Another form of the Schwarzschild metric which is regular 
at the horizon is the Painleve-Gullstrand (PG) form: 

ds^ = - (1 - 2m/r)dr^ - 2^y'2m/rdrdr + dr^ + r'^dVt^ 

= - dr^ + (dr - q^2m/rdT^ + rMfi^ (40) 

(where <^ = —1 means incoming and <^ = means outgoing). The PG time coordinate 

is given by the relation dr = dt — ?^^^^dr. The most noteworthy feature of this form 
of the Schwarzschild metric is that the geometry of the spatial r = constant surfaces is 
flat. We choose the coframe 



= dr, = dr- q^2m/rdT, = rdO, = r sin 9 d(^. (41) 

Some of the connection one-form components, namely F^i, F^i, F^2 have the same values 
as in Minkowski space (so the corresponding AF vanish); whereas the others are 



r^o = -|^2m/r3t9\ F% = <;^2m/rH\ T% = q^2m/rH^ (42) 
(with the corresponding AF having the same values). 
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Epg(N) = -qV2rnrN\ (43) 



For the reference timelike Killing choice, N = dr = Gq — ^y2m/rei, this expression 
yields the value 

Epg{9t) = 2m, everywhere. (44) 

Whereas it is obvious that if N = Bq, then 

Epci^o) = everywhere. (45) 

This latter quasi-local value is consistent with the well-known fact that the ADM energy 
vanishes for the PG metric (since the spatial metric of the constant r surfaces is just 
that of flat Euclidean space). Once again we find the curious result: 

Epcdidr + eo)) =m. (46) 



,2 

On the other hand, for the unit dual mean curvature vector outside the horizon 

1 



N = (eo - c;^2m/rei), (47) 

a/1 — 2m /r 

the PG quasi-local energy has the value 

i^PG(N^) = , , (48) 

W 1 — 2m /r 

with the asymptotic and horizon limits 

Epg(N^),^oo = 2m, ^PG(N^)r=2m ^ oo. (49) 

Note that this value does not approach the ADM energy at spatial infinity. It is well 
known that that desirable property can only be expected to hold for metrics which fall 
off faster than 0(r~^/^), see, e.g., [T5] . 

3.2. FLRW cosmology 

Now let us consider dynamic spherically symmetric metrics. The homogeneous-isotropic 
FLRW cosmological metric has several equivalent manifestly-isotropic-about-one-point 
forms with ds^ = — dt^ + a{t)'^dP, where 

df = dx' + ^\x)dn' (50) 

= (1 -fcr2)-Mr2 + rW (51) 
= 'l + kRy4\~\dR^ + R^dn^) (52) 

= 1 + k{x^ + y^ + z'^)/4\ (da;^ + dy^ + dz^) . (53) 

The first uses the proper radial coordinate p = Xi with S(x) = {sinhx, x, sinx} 
respectively corresponding to the spatial curvature signature k = { — 1,0,+1}. 
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Here we will take the reference metric and connection components to be obtained 
analytically from the respective dynamical ones by taking a{t) = 1, A; = (S = x)? and 
will use the general quasi-local energy expression derived in the Appendix. 

For the first metric form (!50|) . for the quasi- local energy of a sphere at constant t, p 
from f lA15ll with A = a{t), B = a(t)E(x), we find 

^FLRW = -aSAS', (54) 

which is, respectively, 

-Eflrw = a{sinhx(l - cosh x), 0, sin x(l - cosx)}- (55) 

For the area coordinate p = r, from the metric form (15T1) . A = a(l — /cr^)"^/^, B = ar. 
Then from (1A.15P for the quasi-local energy of the sphere at constant t, r we find 



-^FLRW = —arAvl — kr^ = ar(l — vl — kr^). (56) 

For isotropic spherical coordinates take p = R, and from the metric form (152!) 
A = a/[l + {k/4)R'^], B = AR. From (IA.15|) for the quasi-local energy of a sphere 
at constant t,R we find 

akR^ 

^™ = 2[l + ik/A)RY ^^^^ 
We note that the isotropic Cariesmn formula (1A.17P with $ = a[l + {k/4)R'^]~^ obtained 
from the metric form fl53l) gives exactly this same value. 

Although the above results may appear to be different they are in fact identical, as 
can readily be verified using S(x) = r = R/ {1 + kR^/4) with due consideration to the 
respective ranges of the radial coordinates used in these various representations of the 
FLRW metric. In summary, for FLRW we have the respective equivalent quasi-local 
energy values: 



-aR 



3 



Ek=-i = asinhx(l - coshx) = ar[l - Vl + r^] = 2(1 _ R^/^y ' 

Ek=o =0, (59) 

aR^ 



Ek=+i = asinx(l - cosx) = ar[l - - r"^] = 2(1 ^ j^2/4-)2 ' (^0) 

It is noteworthy that, according to this measure, the sign of the quasi- local energy 
is proportional to k, being negative for the open universe, vanishing for the flat case 
and positive for the closed case — but (just as it should be) vanishing when the whole 
universe is considered. 

These results (which were first reported in [16]) may be compared with 
those obtained using the same quasi-local Hamiltonian boundary expression applied 
to homogenous but generally non-isotropic Bianchi cosmological models, using a 
homogeneous choice of reference [17]. That analysis found a vanishing quasi-local value 
for all Bianchi class A models (which includes as special cases the isoropic FLRW k = 
and k = +1 models) and a negative quasi-local energy for all class B models (including 
as a special case the isotropic FLRW k = —1 model). 
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It is also noteworthy that in the FLRW k = —1 model with vanishing matter, one 
finds a{t) = t. It can be directly verified that the geometry is then really fiat Minkowski 
space, yet our quasi-local Hamiltonian boundary term expression gives a non- vanishing 
energy, which, moreover is negative. That a negative quasi-local energy value for certain 
cosmological models can be physically appropriate has been discussed in the work cited 
in the previous paragraph. In the present case, the negative quasi-local energy is related 
to the choice of dynamically expanding comoving observers and their associated choice of 
reference. The next section describes an alternative technique for choosing the reference 
that will yield a different value for the FLRW quasi-local energy in general and for this 
curious special case in particular. 

Remark 

The analytic choice of reference coframe and connection was obtained by taking trivial 
values for some specific constant parameters of the dynamic spacetime. For example, 
for the Schwarzschild case taking m = 0. For this kind of reference choice the quasi- 
local energy may depend on the coordinate systems along the same displacement vector 
(e.g., -^3(1^ ) 7^ Eef(J^ ))• This can happen because this kind of reference choice may 
lead to different reference connections. For example, let the standard Schwarzschild 
reference connection be denoted by Fs and the Eddington-Finkelstein one by Fef- 
It is clear that when taking m = for the reference connection components, both 
of them remain anti-symmetric. If the dynamic orthonormal coframes are related by 
i^EF = Ai^S; then for the reference connection coefficients expressed in the dynamical 
frames Fef ^ -dAA^^ + AFgA-^ 

4. Choice of reference: extremization of energy 

An alternative strategy for obtaining the reference and displacement vector is via 
extremization of the quasi-local energy. We note that Wang-Yau have used this 
technique to select a reference for their quasi-local energy expression [TJ [8]. This is 
reasonable in light of the usual desiderata that quasi-local energy should be non-negative 
and should vanish iff the dynamical space is actually flat Minkowski space. For if one 
supposes that the quasi-local energy expression for any reasonable choice of reference 
indeed were non-negative, and vanished only if the dynamical variables were actually 
those of Minkowski space, then the quasi-local energy could be expected to have a unique 
minimum for some reference. Extremizing the energy w.r.t the reference can be viewed 
as selecting a Minkowski reference that is closest" to the dynamical space (where the 
energy value is used to measure how close). Here we shall assume that the reference 
space has a Minkowski metric and a connection (which, however, need not necessarily 
turn out to be flat; technically we will simply keep the anti- symmetric shape of the 
reference connection when it is expressed in the dynamic orthonormal coframe). 

Since we consider here only spherical symmetric physical spacetimes, the natural 
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choice of the quasi-local two-surface is a constant t, r two-sphere; the tangent space of 
this surface is expressed by the spherical orthonormal frame basis 62, e-^. This simplifies 
the boundary expression and also makes it easier to determine an isometric embedding of 
the two-boundary into the reference space. The extremization comes from extremizing 
the value of the quasi- local energy over the reference gauge choices, i.e., the extremal 
value of energy over the reference coordinate transformations. Through the extremizing 
process, the choice of reference variables and the displacement vector are tied together 
with the dynamic connection. The displacement vector comes out to be the dual mean 
curvature vector of the two-boundary. Using this approach, we are able to obtain a 
quasi-local energy value for the Schwarzschild metric which is independent of the choice 
of the t, r coordinates and which, moreover, gives zero energy for the a = t, k = —1 
FLRW cosmology, i.e., for the dynamic representation of Minkowski space. 

Let us now introduce this process. Suppose the reference metric has the form 

ds^ = -dT^ + di?2 + R^{de^ + sin^ ed$2). (61) 

To determine the quasi-local energy, we have to obtain the reference connection 
which is pulled back from the reference space to the dynamic space via a coordinate 
transformation. This means finding {T, i?, G,$}, which are in general functions of 
{t,r,6,ip}. Because of the special simplicity of the spherically symmetric metrics, we 
can assume the coordinate transformation to have the restricted form 

T = T(t,r), R = R{t,r), = 6, $ = ^. (62) 

Then ( I6T1) becomes 

ds^ = - (T^ - R^)dt^ + 2{RR' - tT')dtdr + (i?'^ - T'^)dr^ + i^Mfi' 

= gQodt^ + 2(7oidtdr + gudr'^ + g22dO'^ + ^733d</^^ (63) 

where t = dT/dt, T' = dT/dr, R = dR/dt, R' = dR/dr. We can rewrite §^ in the 
ADM form: 



^2 



2 



df = -9o^-9^df+ Vmidr + ^dt] +-g22d9' + -gssdv'. (64) 



9ii \ V9u 

Choose the coframe from (16^ 



qo \j9oi- 900911 ^1 ^ goi 

V = y== dt, V = vfi'iidr H — i=dt, 

V9ii V9ii 

32 



= y/^^de = Rde, = Vlasdv? = i? sin ^ dip, (65) 

and the corresponding orthonormal frame is denoted by e^j. Next define the connection 
of the reference: 





= T\ = 


not needed below. 


T\ 


= T\ = 


Pde, T^3 = T\ = P sine dip, 


t\ 


= -t\ 


= Qd9, r^3 = — r'^i = Q sin 9 dip, 


T\ 




= — cos 9 dip. 



(66) 
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where the exphcit form of the component functions are 



P 



Q 



R' 



9oi 



9oo9ii 
-R' 



(67) 



(6^ 



/9ii VR'^ - r'2 

The dynamic connection coefficients have a similar form. The ones that we will 
explicitly need can likewise be parameterized by two functions: 



r"3 



T'^o = Psin^dv?, 



(69) 



1*^2 = r^o = Pd9, 

T\ = = Qde, T\ = -T^i = Qsme dep. 

We assume that the displacement vector = N^Bq + A^^ei is in the normal plane 
of the constant t, r surface, where N^, are functions of {t, r) only, independent of T' 
and R'. The second term of is not vanishing in general. Considering the spherical 
symmetric case, the boundary integral over the constant (t, r) surface S involves the 
Ar^oi term: 



^ j>{D^N^ - D°N^){g22 - 922) sinOde A dip 
lpiArO-DW)((722 



5-22 J 



(70) 



where N^^ is the component expressed in the holonomic basis of reference: N 
iV°eo + A^^ei = N^dr + N^dn. 

The quasi-local energy works out to be 



E 




(71) 



In the previous section, the choices of reference variables for Schwarzschild were 
obtained by taking m = 0. The term / in ( 17T|) vanishes because ((?22 ~ ^22) U = 0. In 
this section, we will let the functions T{t,r) and R{t,r) be undetermined, and through 
the extremization of energy find out what these functions should be. However we will 
require R(t, r)|^ = ^/g22■ This makes the two-surfaces in the dynamic spacetime and the 
reference space isometric. Since the Hamiltonian boundary term is a quantity dependent 
on the quasi-local two-surface, the isometric requirement of the two-surface is reasonable, 
and furthermore, it simplifies the boundary expression by making 1 = 0. 

For any given fixed N, extremize the energy by requiring the vanishing of the partial 
derivative with respect to T' (it is easy to check that taking the partial derivative of 
nl) w.r.t. R' gives the same condition): 



dE 
W' 



R'iT'N^ ± R'N^] 

(_R'2 _ 2^/2^3/2 



0, 



T' 



(72) 
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Substitute into ( ITT]) then 

E = ^[{Q + N^)N'> - {P + N')N']. (73) 

Suppose we choose the normahzed time-hke displacement N which means — (A^°)^ + 
(iVi)2 = then ([73D becomes 

E = v^(l + QN^ - PN^). (74) 



This result implies that the quasi-local energy depends on the free choice of N. We 
further look at the extremal value w.r.t all the displacements. Let A^*^ = cosh a, 
= sinha, take the extremization of the quasi- local energy value (l74l) . 

dE ^ ,^ . , , X n sinha P 

^ = v^(Qsmha-Pcosha)=0, ^ ^ = = (75) 

then we have the relation 

r = T^R', (76) 



and 



Consequently, the quasi-local energy ( ITTl) has the extreme value 



iV° = , ^ = , ^ (77) 

'l - PVQ' Qjl - pyQ^ 



Ee. = v^(l + - ^VQ')- (78) 

We can see that P, Q are determined purely by the metric of the dynamic spacetime. 
There is no longer any information of the reference frame or the displacement vector in 
this energy expression. With the vectors {62,63} tangent to the two-surface, the dual 
mean curvature vector of the two-surface in the dynamic space is N"*" = —keo +pei, 
where k = 2Qj ^fg^ is the extrinsic curvature w.r.t. the space-like normal ei and 
p = —IPjJ'g^ is the extrinsic curvature w.r.t. the time-like normal 69. Then the 



vector N := N-^/|N^| has the components ([77]), where |N-^| := -y/-(N-^, N"^). 
Rewrite (|78l) by replacing Q and P with the extrinsic curvature k and p: 

E.^ = ^ (2/v/^+(-A:)v'1-pV^2) . (79) 

Here we use the definition of extrinsic curvature which is kab '■= —{VaGi.eb), and 
Pab '■= — (Vaeo,e(,), a,b = 2,3; the trace is k = S^-^kab and p = d^-^Pab- By this 
convention k is negative and N"*" is time-like for the dual mean curvature vector, so 



that {—k)yl —p^/k"^ = y/k'^ — p'^ = |N |. Equation fl69|) implies that the trace of the 
reference extrinsic curvature ko is given by taking Q = —1, so that ko = —Ij^fg^. 
Consequently, 

E,. = ^(|N^|-fco), (80) 



2 

which is the same as the Liu-Yau result [9]. Now let us check the following cases. 
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4-1- Standard Schwarzs child. 

The functions necessary here are found from ( fT0l[T2l) to be 

^ = r, P = 0, Q = -^l-2m/r. (81) 

With these expressions, the extreme energy (ITHjl works out to have the standard Brown- 
York value: 



^exs =r(l- v'l-2m/r). (82) 

In this case, as P = we have T' = and the displacement vector N is equal to eg. 
The reference here could be found from T = T{t), R = R{t, r), Q = 6, ^ = (f, with the 
restriction R{t,r)\^ = r. 

4-2. Eddington-Finkelstein. 

The necessary functions obtained from (l29l) now have the values 



, 2m /r ^1 , , 

^=r, P = q ^ ' Q = - (83) 

1 + 2m/ r V 1 + 2m/r 



With these expressions the extreme energy flTHl) again comes out to be 

Ee.EF = r(l - v'l-2m/r), (84) 

which is again the standard value. The condition which restricts the choice of reference 
is 

P 2m 

r = -^R' = c;—R'. (85) 
Q r 

Then one can set any function R[t, r) with the restriction R\g = r, and solve for T{t, r). 
The displacement vector can also be determined from fl77|) to be 

N = (eo - c^— ei)/y'l - 4mVr2, (86) 

which is the dual mean curvature vector. 



4-3. Painleve-Gullstrand. 

The necessary functions found from fl40l) are now 



V^=r, P = q^2mlr, Q = -1. (87) 
Using these expressions the extreme energy fl78p works out to be 



^exPG = r(l- v/l-2m/r). 



Thus once again we found it to have the standard value. The condition which restricts 
the choice of reference is 
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The displacement vector is N = (eo — <;^2m/rei) / ^1 — 2m/r. 

In a similar fashion, if one considers the functions ^/g^, P, and Q associated with 
the spherical isotropic Schwarzschild coframe fl22l) . one will once again obtain from fl78|) 
the standard quasi-local energy value. 

4. 4- FLRW cosmology 

For the dynamic FLRW cosmological models, from the extreme energy expressions flTHl) 
it is sufficient to calculate the quasi-local energy using the metric functions ^Jg^ = a{t)r, 



P = ar, Q = —y/l — kr"^ obtained from (!50|) — the other forms of the FLRW metric 
would lead to the same answer. The result is 



ar'^[k + a^^ 



i^exFLRW = ar(l - Vl - kr^ - a^r^) = . ' , ^ ^ ^ . (90) 

In contrast to the analytic result (156|) . which can be negative, with the aid of the 
Friedmann cosmological equation, 
k Stt 

it can be seen that this value is non-negative. There is no contradiction, the present 
quasi- local energy value (pOj) corresponds to non- expanding observers using a different 
reference. 

Consider in particular the special test case k = —1, a{t) = t. We find the energy 
value: E' = 0. We can find a simple reference choice R = tr (actually, this needs be 
satisfied only on the two-boundary, not necessary in the whole space), then from (1761) 
we obtain T = tA/l + r^. It is well-known that the 4-geometry 

- dt^ + ^dr^ + t^rMn'^ (92) 

1 -\- r'^ 

is actually Minkowski space, and zero energy is just the value we expect. 



4-5. Discussion 

In contrast to the analytic approach, here we used energy extremization to select the 
reference and displacement vector field. We tested the resulting quasi-local energy 
expression on several forms of the Schwarzschild metric, obtaining in each case the 
standard quasi-local energy value. We also tested the expression on the FLRW 
cosmological metric, obtaining a new result for the FLRW quasi-local energy. In 
both cases the time displacement vector field turns out to be the dual mean curvature 
vector, and the quasi-local energy value has the desirable property of being non-negative, 
vanishing iff the dynamic geometry is fiat Minkowski space. 



5. Conclusion 



The covariant Hamiltonian formalism has been incomplete in one aspect. The 
Hamiltonian boundary term, whose value determined the quasi-local quantities, in 
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addition to depending on the dynamical fields, also necessarily depends on reference 
values for these dynamical fields (which specify the ground state with vanishing quasi- 
local quantities) along with a spacetime displacement vector field. However, no specific 
proposal had been made as to how to choose these unspecified quantities. Here, for 
Einstein's GR, for certain spherically metrics (the static Schwarzschild metric and the 
dynamical homogeneous isotropic cosmologies), following [10] we have considered two 
techniques for choosing the reference and vector field. 

The first (which goes back to [Hill]) depends on an analytic choice of reference fields 
obtained by taking trivial values for certain parameters in the metric and connection 
coefficients. For the usual time slicing for several spatial metrics this leads to the 
standard Brown- York quasi-local energy for the Schwarzschild metric, but to different 
energy expressions for the alternative time slicings of the Eddington-Finkelstein and 
Painlave-Gulstrand metrics. For the FLRW cosmological metrics it leads to a quasi- 
local energy value which is proportional to the sign of the spatial curvature (and thus a 
negative energy for a certain dynamical slicing of Minkowski space). 

The other approach uses extremization of the quasi-local energy to select an optimal 
reference and timelike vector. The resulting quasi-local energy for these spherically 
symmetric metrics is independent of the coordinates and is non-negative (for both the 
Schwarszchild and FLRW metrics) and vanishes only for Minkowski space. For the 
Schwarzchild metric it gives the standard quasi-local value. Going beyond the present 
work, there have been further developments in the energy optimization approach; a brief 
letter describing this has already appeared [18] . 

Appendix A. General formulas for spherical analytical reference choice 

Here we briefly present the quasi-local energy calculations using an analytic reference 
choice for general spherically symmetric metrics in several different coordinates and 
coframes. 

Spherical frames 

Consider orthonormal co-frames using spherical coordinates of the form: 



where N, A, B are functions only of the general time and radial coordinates r, p. The 
(metric-compatible hence anti-symmetric) connection one-form coefficients are readily 
obtained from the differentials 



■=N dr, := A dp, := BdO, d'^ := B sin 9 df/?. 



(A.l) 



d^p 



{N'/NA)'&P 
(A/AN)^'' A i^P, 



(A2) 
(A.3) 



d& 



{B/BN)d^ A + {B'/AB)dP A 



(A.4) 



di9^ 



{B/BN)§^ Ad"^ + {B'/AB)r A + (1/5) cot 9^^ A ^9^, 



(A.5) 
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where dot and prime represent respectively the r and p partial derivatives. The ones of 
particular interest to us are 

r% = -^^^ = ^de, r% = = f ^in ^ d^. (a.6) 

We will also need the associated reference values, 

f % = (B'/A) de, r% = (B'/A) sin e d^, (a.t) 

which we here have assumed to be given analytically (by taking limits like m — )■ 0) and 
thereby affecting the transformations A A, B ^ B. In our calculation we will need 

Ar% = A{B'/A) de, AT^p = A{B'/A) sin 6 d^. (A.8) 
Cartesian frame and coordinates 

Spherically symmetric metrics may be also be described by Cartesian coframes using 
Cartesian spatial coordinates (labeled by latin indies with range 1,2,3) in the form 

^^ = NdT, ^9'^ = $dx^ (A.9) 

where A^, $ are functions of t, R with = x^Xk- We find 

d^9" = {N'/N^R)xk^^ A (A. 10) 

d^'' = {^/N<t>)^^ A^'' + /^'^R)xm^'^ A d^. (A.ll) 

The connection coefficients of particular interest are found to be 

pi _ ($7<|.2i^)(2;%^ _ x'^^). (A. 12) 

We here assume that the associated reference is given by the Minkowski space obtained 
by analytically lestiicting these formulas to = 1, $ = 1; thus the reference connection 
is such that its values vanish. This is just what we expect for the Minkowski frame 
determined by the coordinates x'^. Since the reference connection vanishes the Cartesian 
case provides a good check for our calculations. 

Our choice of particular non-vanishing reference values for the various spherical 
representations can be understood as just what is needed, as we shall see, to arrange to 
give the same results in all of these frames. 



Energy expression 

We are interested in the particular quasi-local energy given by our preferred Hamiltonian 
boundary term 2-form expression ([9]) with vanishing 2nd term: 

2kB{'N) := Ar°/3 A L^7]J. (A.13) 

Here we will take N to be the unit time-like displacement on the boundary, which is at 
constant p, r. The other choices considered in the text are proportional to this choice 
For spherical frames with the displacement choice N = e_|_ (i.e., one unit of proper 
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time orthogonal to the constant "time" hypersurface), our Hamiltonian boundary term 
quasi-local energy 2-form expression reduces to 

2KBie^) = AT'"' A r]^ab = A v^^e + 2Ar^^ A 7]^^^ 

= AAT^^ A r]^pe = 4Ar^^ A^"^ = -4BA{B'/A)dQ. (A.14) 

The associated quasi-local energy, obtained by integrating over a 2-sphere at constant 
r, p (with K = Stt), has the value 

Es{e^) = -BA{B'/A). (A.15) 

It is notable that the result is not explicitly dependent on the lapse N. 
On the other hand, for the Cartesian frame we find 

2KB{e^) = Ar^ A 7]^ij = 2($7$2i?)a;i^» a r/^ij 

= - 4($7$2i?)a;Vfc = -4*'i?^dl]. (A.16) 

Hence, in this case, the quasi-local energy obtained from integration over the 2-sphere 
at constant r, R is given simply by 

^c(e±) = (A.17) 
Schwarzschild Applications 

In particular we have for the area coordinate from f lTT]) A = {I — 2m/r)~^/^, B = r. 
Using these in the general spherical energy expression ( 1A.15I) gives 

^s(e±) =r[l-(l-2m/r)^/2], (A.18) 

which is the standard energy value ( IT8|) . An equivalent expression, 

Esie^) = . (A.19) 

1 + Jl- 2m/r 



more clearly reveals the horizon and asymptotic limits. For isotropic spherical 
coordinates, from (El]), A = {I + m/2Kf, B = R{1 + m/2Rf. Using these in the 
general spherical expression flA.lSP yields the quasi-local energy 

Eg(e^) = m(l + m/2R). (A.20) 

Recalling that r = R{1 + m/2R)^, we find that this is actually the same as the standard 
value f lA.18|) . On the other hand for the isotropic Cartesian frame, using from ( 128|) . 
$ = (1 + m/2RY in the Cartesian expression (IA.17p . turns out to give the same value 
(IA.20p . equivalent to the standard value (lA.lSp . As this case has a vanishing reference 



connection, it provides an important confirmation, not only for the standard energy 
value, but also especially for our analytic technique of choosing the the reference. 
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